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The various flow phenomena observed at energies between 50 AMeV and 160 AGeV 
are reviewed. I first define three types of flow: directed flow and elliptic flow, which are 
the two first Fourier components of the azimuthal distribution in non-central collisions; 
radial flow, which is deduced from an analysis of transverse momentum spectra in central 
collisions. Then, I review the observations of directed flow and elliptic flow, with emphasis 
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on recent results. I discuss their dependence on various parameters: global geometry (im- 
pact parameter, mass numbers of colliding nuclei and bombarding energy) and individual 
observables (rapidity, transverse momentum and particle type). Finally, I explain how 
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azimuthal distributions can be measured experimentally. 



^ : 1. THREE TYPES OF FLOW 

The azimuthal angles of particles emitted in a nuclear collision are correlated with 
the direction of impact parameter In this talk, I use the generic term "flow" to 
refer to these correlations (which are essentially macroscopic effects), irrespective of any 
theoretical interpretation. 

The identification of the direction of impact parameter in a high energy collision is a 
non trivial task, which relies on a study of azimuthal correlations between the reaction 
products (see Section 4). In practice, this identification is possible only if the system is 
large enough: in proton-proton and proton-deuteron collisions at 28 GeV, the observed 
azimuthal correlations are back to back and can be explained by momentum conservation 
alone @. In proton-nucleus collisions, flow may be present, as indicated by a recent study 
of correlations between pions and nucleons ||. However, flow is seen unambiguously only 
in nucleus-nucleus collisions, where it was discovered at Bevalac in 1984 |l|. It is thought 
to be highly sensitive to medium properties, and has therefore been extensively studied 
in recent years. 

In this section, I denote by the azimuthal angle of an outgoing particle, measured from 
the direction of impact parameter (see Fig.|I]). There is flow if the distribution is not 
isotropic. A convenient measure of flow is provided by the Fourier coefficients, (cosn</>) 
(where the brackets denote an average value in a given kinematic {pt,v) window). The 
two major flow effects observed so far correspond to the n = 1 and n = 2 | Fourier 
components. These effects, which I shall name directed and elliptic flow respectively, are 
illustrated in Fig.|l]. In the projectile rapidity region, directed flow is positive if (cos0) > 0, 
and negative if (cos0) < 0. For symmetric collisions, (cos0) is an odd function of the 
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Figure 1. The four major types of azimuthal anisotropies, viewed in the transverse plane. 
The target is denoted by T and the projectile by P. Top: Directed flow in the projectile 
rapidity region, positive (left) and negative (right). In the target rapidity region, the left 
and right figures are interchanged. Bottom: elliptic flow, in-plane (left) and out-of-plane 
(right). 



centre-of-mass rapidity, and signs are therefore reversed in the target rapidity region. 
Similary, for elliptic flow, I distinguish in-plane elliptic flow if (cos 20) > 0, and out-of- 
plane elliptic flow if (cos 20) < 0. Note that unlike directed flow, elliptic flow has the 
same sign in the target and projectile rapidity regions, at least for symmetric systems. 

In central collisions, the direction of impact parameter cannot be defined: azimuthal 
distributions become flat, and directed and elliptic flow vanish. However, a third type 
of flow, called radial How, has been defined for central collisions (although it is not a 
"flow" in the sense defined above) |5|-|7|. The idea is to compare the measured kinetic 
energy (or transverse mass in relativistic collisions) spectra to those of a thermal fluid 
expanding outwards Consider, for simplicity, a non-relativistic perfect fluid with 

temperature T and velocity Vfl ow . The velocity of a particle in the fluid can be written as 
v = Vfl ow + v t h, where v t h is the velocity of the random thermal motion. Then the average 
kinetic energy of a particle with mass m is (E) = mvl ow /2 + 3T/2: while the thermal 
energy is independent of m, the collective flow contribution is proportional to m. From 
the spectra of particles with different masses, one can thus extract both the thermal and 
the collective component. This will not be discussed further here. 

Fig|J displays the typical variation of directed and elliptic flow with the bombarding 
energy. The high energy part will be discussed in detail in the next two sections. Here, 
I comment on the change of sign of both directed and elliptic flow below 100 AMeV. At 
low energies, the interaction is dominated by the attractive nuclear mean field, which as 
two effects: first, projectile nucleons are deflected towards the target, resulting in nega- 
tive directed flow jlOfl ; second, the projectile and target form a rotating system, and the 
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Figure 2. Schematic behaviour of the magni- 
tudes of directed flow (top) and elliptic flow 
(bottom) function of the bombarding ki- 
netic energy per nucleon in the laboratory 
frame. Full lines: proton flow; dashed lines: 
pion flow. 



Figure 3. Energies of vanishing directed 
flow (E B al) and elliptic flow (E TR a), as a 
function of impact parameter (from [15]), 
in the centre of mass system. The corre- 
sponding laboratory energies are larger 
by a factor of 4. 



centrifugal force emits particles in the rotation plane ("rotation-like effect" JTTJ] ) , produc- 
ing in-plane elliptic flow [TI|. At higher energies, individual nucleon-nucleon collisions 
dominate over mean field effects. They produce a positive pressure, which deflects the 
projectile and intermediate rapidity fragments away from the target ("bounce-off" and 
"sidesplash" effects [ft]), resulting in positive directed flow. Furthermore, the participant 
nucleons, which are compressed in the region where the target and the projectile over- 
lap (see Fig.|TJ) , cannot escape in the reaction plane due to the presence of the spectator 
nucleons ("squeeze-out effect" Q), producing out-of-plane elliptic flow. 

The crossing energies at which directed [13] and elliptic [Q flow cancel are displayed in 
Fig. 3 for Au-Au collisions, as a function of impact parameter. An extrapolation to zero 
impact parameter gives values which are consistent, within error bars, with the threshold 
energy for radial flow, estimated as -Erad = 8.7 ± 2.5 AMeV. As Crochet et al. point out 



151 , this suggests that a common mechanism is at the origin of positive directed flow, 
out-of-plane elliptic flow and radial flow: the three phenomena appear when the attractive 
forces are counterbalanced by the thermal pressure. A possible relation of this crossing 



phenomenon with the nuclear liquid-gas phase transition is discussed in |16 . 

I now turn to a more detailed analysis of directed and elliptic flow at energies above 
100 MeV per nucleon. 
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2. DIRECTED FLOW 



The positive directed flow of nucleons has been extensively studied below 2 GeV per 
nucleon, both theoretically and experimentally Since it is essentially a pressure effect, 
it is naturally expected to be sensitive to the compressibility of nuclear matter. The 



momentum dependence of nuclear interactions is also important [17 



The magnitude of directed flow depends on the centrality of the collision: it vanishes 
for central collisions by symmetry, and for very peripheral collisions where compressional 
effects are low. It is therefore maximum for semicentral collisions 
parameter of the order of 5-6 fm for Au-Au collisions 



18|, at an impact 



19 . The measurements of directed 



flow to be discussed below are made at the centrality where it is largest. 

As mentioned in Section 1, directed flow is an odd function of the centre-of-mass rapidity 
in symmetric collisions. It is therefore linear near mid-rapidity. Furthermore, a saturation 
is observed near projectile and target rapidities, resulting in a typical S shape |2(J . 

Due to this peculiar dependence, and to acceptance considerations, the strength of 
directed flow is usually measured by plotting the average momentum projected on the 
direction of impact parameter, (p? cos cj)), as a function of rapidity, and by taking the 
slope of this curve near mid-rapidity, the so-called flow parameter [|T8|j : 
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dy 



(1) 



y=yc 



where y pTO j denotes the projectile rapidity. 

The flow parameter of protons scales with the mass numbers A and B of the colliding 
nuclei approximately like A 1 / 3 + B 1 / 3 



This is an important point. If the system 
behaved like a perfect fluid, it would obey geometric scaling [^2| and the value of F would 
be independent of A for a symmetric collision (neglecting the nuclear skin depth). The 
fact that F is proportional to A 1 ! 3 + B 1 ^ 3 , i.e. to the collision time, which was predicted 
by microscopic calculations El, is an indication that the system is on its way towards 



thermalisation, but only partially thermalised [p4 



Further insight into the mechanism of directed flow is obtained by studying the pro- 
duction of fragments. Their directed flow increases with their mass at Bevalac and SIS 
energies JI^|P7[. A similar behaviour at AGS energies was reported at this conference 



[34]. Repeating the simple fluid dynamical argument given in Section 1 for radial flow, 
this suggests that directed flow is an effect of collective motion in the expanding system. 

The variation of directed flow with bombarding energy is also instructive: between 
150 MeV and 1 GeV per nucleon, i.e. in the non relativistic regime (recall that energies 
in the centre of mass are 4 times smaller), the flow parameter F is approximately propor- 
tional to the projectile momentum ||25|| . In other terms, the velocities of outgoing particles 
scale like the incoming velocity t>j. Although this scaling was predicted in a hydrodynam- 
ical model [p2], it is not an effect of thermalisation. It simply reflects the fact that in 



this energy range, the only relevant velocity scale is effects associated with the Fermi 
velocity Vp or with relativity are negligible in the energy regime where v p *C <C c. 

At lower energies, attractive effects become important, and the scaling breaks down. 
Furthermore, since the attractive part is essentially a surface effect, the A dependence 
becomes more complicated f26fl . 
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At relativistic energies, the behaviour is also different. The scaled flow F s = Fj (A 1 / 3 + 
B 1 / 3 ), whose value below 1 AGeV is approximately F s ~ 0.03p pro j, reaches a maximum 
of 50 MeV and then decreases slowly. A first indication of this decrease between 1 and 
2 AGeV was reported at the last Quark Matter conference |[25||. It is confirmed between 



2 and 6 AGeV by results presented at this conference ||27| . At 10.8 AGeV, the E877 
collaboration measures F s ~ 35 MeV |28||, compatible with the value F s ~ 31 ±7±6 MeV 



found by E917 [29 



The first evidence of directed flow at SPS was reported by the WA98 collaboration 
at the Jaipur Conference this year [30|. More recent measurements have been made by 



the NA49 ^ and NA45 H collaborations. Directed flow seems to be concentrated m 
the fragmentation regions, and almost absent in the central rapidity region: the rapidity 
dependence does not follow an S-shaped curve (deviations from this shape are already 
seen at AGS by E917 |^|) and the flow parameter, defined as the slope near mid-rapidity 
(Eq.(|lD), becomes irrelevant. The value of (prcos0) m the fragmentation region is at 
least three times smaller than at AGS |pC|| , confirming the decrease discussed above. 

The flow parameter F is a pT~integrated observable. A more detailed information on 
directed flow can be obtained by plotting (cos0) as a function of p?- At low p T , (cos</>) 
varies linearly with p^- This is because the momentum distribution must be a regular 
function of p x = pr cos <p and p y = px sin (p. It was recently pointed out that 
directed flow can be negative at low pt if strong radial flow occurs in the system. At 
higher pt, the variation of directed flow gives an information on momentum dependent 
interactions |35|,|36| . 

Directed flow of particles other than nucleons has been observed, but the available 
information is less detailed. Pions exhibit weak negative directed flow (i.e. opposite to 
the nucleons), first observed in asymmetric collisions |38|j39[1 . This is interpreted as a 
consequence of pion rescattering on spectator matter (ttN — > itN) P0 |. The same effect 
has been recently observed at AGS j|T| (although the flow is positive for high p?, high 
rapidity pions) and at SPS |3|, |30|j3T|| . 

Directed flow of strange particles has been also studied at Bevalac and SIS. While A 



baryons exhibit positive directed flow, of roughly the same magnitude as protons [42 



kaon directed flow is compatible with zero j|3J. This may seem surprising since only 
associate production is kinematically allowed at these energies (2 AGeV), and A's and 
fT's are therefore produced in pairs. Since K +, s have a large mean free path in nuclear 
matter (unlike pions, they cannot be absorbed), their flow provides a clean probe for the 
kaon potential in a nuclear medium HPfl - First data on K + directed flow at AGS were 



presented at this conference [34]. The value is also consistent with 0, except for low-pr 



kaons where it is negative, i.e. opposite to the proton flow, which can be interpreted as 
an effect of Coulomb repulsion. 

The first observation of antiproton directed flow was reported at this conference [45], 
but the sign is unclear yet. Negative flow is expected, due to pp annihilation 

Finally, the isospin dependence of directed flow may give interesting results: Coulomb 
effects, among others [[|7Jj, are expected to produce a small isospin dependence which, if 
measured, could provide a direct information on the size of the system. Experimentally, 
this can be studied in two different ways: first, one can compare the magnitude of proton 
directed flow in colliding systems with different isospins. This has been done recently 
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at low energies 1 4~8fl. Alternatively, one can compare the magnitude of directed flow for 
particles of a given isospin multiplet. No obvious difference is seen between neutron and 
proton [|{J, or between K + and K® 
low transverse momenta is seen by E877 BT 



4^]. A difference between ir + and it flows at very 



3. ELLIPTIC FLOW 



While directed flow keeps the same sign from SIS to SPS energies, a transition from 
out-of-plane to in-plane elliptic flow occurs at relativistic energies, which was predicted a 
few years ago |50|-p^] and recently observed pHf. In this section, I first recall a few general 



properties of elliptic flow. Then, I review the observations at SIS energies, and the recent 
measurements at AGS and SPS. 

The impact parameter dependence of elliptic flow is similar to that of directed flow [|55 



However, the maximum is found at a slightly higher impact parameter, about 7 fm for 
j54|, instead of 5-6 fm for directed flow. On the other hand, contrary to 

At SPS, 



Au-Au collisions 



directed flow, the rapidity dependence is essentially flat at SIS and AGS |p5jp8 



a maximum is observed at mid-rapidity for protons pT |. 

Elliptic flow is higher for fragments than for protons @-|56|: as for directed flow, this 
is a signature of collective motion. However, the magnitude of elliptic flow is larger for 
larger nuclei |55| , [57|| , which indicates that only partial thermalisation is achieved. 

Elliptic flow is always out-of-plane at Bevalac and SIS energies. It reaches a maximum 
at about 400 AMeV p5|,|56|,|58| . However, the variation between 250 and 800 AMeV is 
very weak [51], which suggests the same scaling behaviour as observed for directed flow in 
this energy region. This scaling can also be seen by studying the pt dependence of elliptic 
flow. At low pt, simple mathematical arguments (similar to those used for directed flow) 
show that it is proportional to p\ |59]]. Then it rises as a function of pt, and a universal 
curve is obtained when p? is scaled by the projectile momentum p W o] |54,58|. 

Charged |JIJ and neutral [5T| pions also undergo out-of-plane elliptic flow. The dom- 
inant contribution to this effect is absorption on nucleons [ttNN — > AN — > NN) [|62jl . 
Concerning the isospin dependence, the magnitudes are similar for neutrons and protons 

MM 



and no difference is seen between n + and 7T~ 
The transition from out-of-plane to in-plane elliptic flow at relativistic energies is a 
consequence of the Lorentz contraction of the spectators. They leave the interaction 
region after a time of the order 2R/ r y (R being the nuclear radius and 7 the Lorentz 
contraction factor). Later, outgoing particles are free to escape anywhere in the transverse 
plane. However, if partial thermalisation occurs, pressure will drive them along the lines 
of pressure gradient (recall that the force per unit volume is —VP). Since the overlap 
region between target and projectile has a smaller size along the impact parameter (see 
Fig.l), the pressure gradient is larger in the reaction plane ||50|| . This produces in-plane 
elliptic flow. 

First indications of non-zero elliptic flow were obtained a few years ago at AGS by E877 
[EI| and at SPS by WA93 and NA49 However, conclusive evidence that it was 



in-plane came only recently | 28fl . At AGS, in-plane elliptic flow is seen for protons, but 



not for pions j|T|, while it is seen for both pions and protons at SPS |30| , |3T 



From the qualitative arguments given above, one expects that the transition to in- 
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plane elliptic flow occurs when the Lorentz contraction becomes significant ||5T| . The 
first measurement of the transition energy was reported at this conference by the E895 



collaboration, who obtain a value lying between 4 and 6 AGeV f27j| . At these energies, 
the kinetic energy is almost equal to the mass energy in the centre of mass system, thus 
confirming the validity the above argument. A more thorough analysis has recently shown 
that the transition energy reflects the time dependence of the pressure [p8|J . 

Let us finally comment on the magnitude of in-plane elliptic flow: the integrated value 



for protons increases approximately from 2% to 6% between AGS [28] and SPS [31 



RQMD calculations |69| predict values of the order of 4% at SPS, in qualitative agree- 



ment with experimental data |31|j32][ , but significantly lower than the predictions of a 
hydrodynamical calculation pjQfl . This suggests that that the colliding system is not fully 



thermalised locally. The variation of in-plane elliptic flow with the number of collisions 
per particle, i.e. with the degree of thermalisation achieved in the system, is studied in 

roll. 



4. FLOW ANALYSIS 

I now briefly explain how the results presented in the previous sections can be obtained 
experimentally. First, one must estimate the direction of impact parameter event by 
event. If positive directed flow is present, adding the transverse momenta of all particles 
detected in the projectile rapidity region yields a vector which points in the direction of 
impact parameter (see FigJ]]). This forms the basis of the transverse momentum method 
ZD|. More generally, the azimuthal angle of the impact parameter (measured from 



a fixed point in the detector) is estimated by constructing the following vector in the 
transverse plane: 



Q = ( n C 2tl ) = X>* ( C °*t ) P> 



Qsm<f> R J ^ K \sm 



>k 



where the sum runs over detected particles with azimuthal angles <pk, and Wk is a weight 
which is positive (resp. negative) in the projectile (resp. target) rapidity region. Wk can 
be any function of pr and y, and must be optimised to yield the best possible reaction 
plane determination. 

The direction of impact parameter, estimated from Eq.(Q), differs from the true one by 
an error A0, due to finite multiplicity fluctuations (see Fig.^). The effect of this error on 



flow measurements must be corrected. I now explain a method [|7l|] which allows to do 
this correction systematically. 

The measured azimuthal angle ip of an outgoing particle is related to the true azimuthal 
angle <fi by ip = — A<f>. Averaging over many events, assuming that <fi and A0 are 
independent (autocorrelations are avoided by removing the particle from the sum in Eq.(^j) 
70[), one obtains the following relation between measured and true Fourier coefficients 



urns 



(costj/0) = (cosn0)(cosnA0). (3) 

From Eq.(|3p, one can reconstruct the true distribution once the correction factor (cos nA<f)) 
is known. 
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Figure 4. Schematic picture of the distri- Figure 5. Variation of (cosnA0) with 

bution of Q, given by Eq.(4). Q fluctuates the parameter x, calculated from Eq. (5). 

around its average value (Q) = Q e x with a The curves are labeled by the value of n. 
standard deviation a, so that a typical value 
of Q (dotted arrow) lies within the shaded cir- 
cle with radius a. 



The important point is that, as we shall see shortly, the distribution of A0 is a univer- 
sal function of a single dimensionless parameter Xi which measures the accuracy of the 
reaction plane determination and can be determined directly from experiment. 

In order to calculate the distribution of A0, we must parametrise the fluctuations of 
Q. Since the number N of particles entering the definition of the vector Q in Eq.@ is 
usually much larger than unity, the central limit theorem can be applied. It shows that, 
for given magnitude and orientation of the true impact parameter, the fluctuations of Q 
around its average value (Q) = Q e x (which is the true direction of impact parameter) 
are gaussian. Then the two dimensional distribution of Q takes the form |50|,|52|,[73|,[f4[ 



dN 1 ( 1Q-(Q)1 2 ^ 1 ( Q 2 + Q 2 -2QQcosA(J) \ 



QdQdA(f) no 2 r V a 2 J na 2 ^ \ a 2 J 



Note that Q scales like N while a scales like v N. 

Eq.(||) can be integrated over Q |52| to yield the distribution of A<j>. As can be seen 



in Fig.|], the typical magnitude of A0 is a/Q. It is therefore useful to introduce the 
dimensionless parameter x = Q/ a : which characterises the accuracy of the reaction plane 
determination. The correction factors (cos?iA0) entering Eq.(H) depend on Q and a only 
through x- They are calculated by integration of Eq.(|4]) over A</> and Q [ [Tl[] : 

(cosnA0) = ^ X e-* 2/2 1^1 (^) + % {^) (5) 

where Ik is the modified Bessel function of order k. The variations of the first coefficients 
with x is displayed in Fig.[5]. Since A0 is of order (cos nA<j)) is close to unity for large 
X- 
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The last step of our reconstruction procedure is to estimate the value of x-, i- e - to 
measure the accuracy of the reaction plane determination. The most popular method 
to do so [19,42 ,|54 , |60 , |61 is to divide each event randomly into two subevents containing 
half of the particles each, and to construct Q independently for each subevent One 
thus obtains two vectors Qj and Qjj. The distributions of Qj and Q /7 are given by an 
equation similar to Eq.(^j). However, since each subevent contains only N/2 particles, 
the corresponding average value and fluctuations must be scaled: Qj = Qjj = Q/2, 
&i = &11 = o~/V2, and therefore xi — Xn — x/ The distribution of the relative angle 
between Q 7 and Q //; A(p R = |A0j — A<p n \ can be calculated analytically [|5^,[74[]. The 
fraction of events for which A<p R > 90° is simply related to %: 



iV(90° < Acj) R < 180°) _ exp(-x 2 /2) 
iV(0° < Acj) R < 180°) ~ 2 ' 



(6) 



Once x is known, the Fourier coefficients of the azimuthal distribution can be recon- 
structed from Eqs.([|) an d @. Until now, only the two first components, n — \ (directed 
flow) and n = 2 (elliptic flow) are measured. It would be interesting to extend these mea- 
surements to higher harmonics [73]. Upper bounds on the next two Fourier coefficients 
have been given by E877 [|IJ,BS|. 

The transverse momentum method, Eq.(|2|), gives a reasonable estimate of the reaction 
plane if directed flow is large enough. However, directed flow is weak at SPS and will 
probably not be observable at RHIC and LHC. For this reason, an alternative method 



to determine the reaction plane was proposed in [52], which makes use of elliptic flow, 
instead of directed flow. If in-plane elliptic flow is present, is estimated according to 
the following equation: 



Q 2 



Q 2 cos2$ R 
Q 2 sin2$fl 



N 



k=l 



COS. 

sin 



(?) 



In this case, the weight is positive everywhere, since the sign of elliptic flow does 
not depend on rapidity. This method is equivalent to a diagonalisation of the transverse 
sphericity tensor ||50|| , and was proposed independently in for low energies where 
directed flow also vanishes. Note that the last equation allows to determine $r modulo 
7i. The magnitude of elliptic flow at SPS is determined more accurately with this method 
than with the transverse momentum method 31 . 

Note that the above methods do not allow to determine the sign of the flow. Eq.([|) 
assumes that directed flow is positive (if it was negative, should be replaced by $r+7t 
in Eq.(|2])). Similarly, Eq.(|7]) assumes that elliptic flow is in-plane. In fact, the sign of 
directed flow was determined only at low energies, where it is negative, by coincident 
polarisation measurements fllQfl . The sign of elliptic flow can be determined only by a 
simultaneous measurement of directed flow. 

One must be aware that the above description of flow analysis is oversimplified. In 
particular, Eqs.@ and (0) cannot be used in experiments where particles are not seen 
individually, but only through the energy deposited in a calorimeter, as in most AGS and 
SPS experiments. I refer the reader to the corresponding publications for more details 
28|,|1],|I],||,|7|. 
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5. CONCLUSIONS AND OUTLOOK 

Directed and elliptic flow are present at all energies in nuclear collisions. Their mea- 
surement provides direct insight into final state interactions, which can only be deduced 
indirectly from other observables such as inclusive distributions. Therefore, flow is highly 
sensitive to medium properties. Detailed experimental analyses, carried out with large 
acceptance detectors, are now available at Bevalac and SIS energies. Since a given ob- 
servable usually gets contributions from several effects, it is difficult in practice to extract 
accurate numbers for physical quantities such as compressibility, optical potentials, in- 
medium cross sections. Nevertheless, it has been recently emphasized that no model is 
presently able to explain all the data quantitatively |19|| , which shows that flow provides 
useful constraints on theory. 

The increase of directed and elliptic flow with fragment mass indicates that they reflect 
the collective motion of the expanding system. On the other hand, their increase with 
the mass numbers of colliding nuclei indicates that there is no geometric scaling at SIS 
energies, which implies that thermalisation is only partial. The small magnitude of elliptic 
flow at SPS suggests that the same holds at SPS. Flow measurements in collisions of 
lighter nuclei could provide some insight into the highly debated issue of thermalisation 
at ultrarelativistic energies. 

A wealth of new results have been obtained since the last Quark Matter conference. 
Flow is now studied at AGS by E877 g|, E895 J27[ and E917 @, and at SPS by NA45 



NA49 [0, NA52 [0 and WA98 |JU|. Directed flow and in-plane elliptic flow are now 



clearly seen at AGS and SPS, and flow of identified particles is measured p7|j29| - pTl , p| . [15 



It would be interesting to extend these measurements to other particles such as multi- 
strange baryons, dileptons, J/^'s. This would further constrain enhancement /suppression 
mechanisms. 

While directed flow becomes smaller as the bombarding energy increases, in-plane el- 
liptic flow becomes larger. One can reasonably expect that the latter will be the only 
observable flow effect at RHIC and LHC. Since in-plane elliptic flow results from com- 
pression of the produced hadronic matter, its magnitude should give information on the 
pressure, i.e. of the equation of state. 

In order to obtain a consistent description of heavy ion collisions, flow analyses should 
eventually be correlated with other observables. For instance, the knowledge of the reac- 
tion plane should allow to obtain a more detailed information on the space-time distri- 



bution from HBT, 2 particle correlations [[76l , [77| . A unified description of radial, directed 



and elliptic flow could also shed some light on the issue of thermalisation p3| 
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